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ABSTRACT 

Observations with the Extreme Ultraviolet Explorer (EUVE) have shown the Coma Cluster 
to be a source of EUV emission in excess of that produced by X-ray gas in the cluster. We have 
re-examined the EUVE data on this cluster in an attempt to obtain clues as to the origin of this 
emission. We find two important new results. First, the ratio between the azimuthally averaged 
EUV excess emission and the ROSAT hard X-ray flux is constant as a function of distance from 
the cluster center outward. Second, a correlation analysis between the EUV excess emission and 
the X-ray emission shows that on a detailed level the EUV excess is spatially closely related to 
the X-ray emission. These findings contradict previous suggestions as to the underlying source 
of the diffuse EUV emission in Coma and provide important information in regards to the true 
source of this emission. We propose a new explanation for the source of this emission: inverse 
Compton scattering of microwave background photons by secondary electrons and positrons. 
We explore this possibility in some detail and show that it is consistent with all of the available 
observational evidence. The parent cosmic ray protons may have been produced by any of a 
number of sources, including supernovae, active galaxies, galactic winds, and cluster formation 
shocks, but we believe that the most likely source is cluster formation shocks. If the EUV 
emission in the Coma Cluster is, in fact, the result of secondary electrons, this may be the only 
direct evidence for secondary electrons in the intracluster medium of a cluster of galaxies, since 
recent work suggests that secondary electrons may not be the cause of radio halos. 



1. Introduction 

Observations with the Extreme Ultraviolet Ex- 
plorer (EUVE) provided evidence that a number 
of clusters of galaxies emit excess EUV emission 
in the cores of the clusters. The first clusters re- 
ported to have EUV excesses were the Virgo clus- 
ter (Lieu et al. 1996a; Bowyer et al. 1996) and 
the Coma Cluster (Lieu et al. 1996b). Thereafter 
EUV emission was reported for Abell 1795 (Mit- 
taz. Lieu & Lockman, 1998) and Abell 2199 (Lieu 
et al. 1999a). These early works employed a vari- 
ety of data analysis schemes that were later found 
to be incorrect (Bowyer, Berghofer, & Korpela, 



1999), primarily because incorrect methods were 
used to account for the sensitivity profile, or expo- 
sure map, of the telescope. The only clusters that 
have been determined to have an EUV excess us- 
ing uncontested data analysis procedures are the 
Virgo cluster (Berghofer et al. 2000) and the Coma 
Cluster (Bowyer et al. 1999). 

Subsequent to the analysis of Bowyer et al. (1999) 
additional EUV data on the Coma Cluster were 
obtained with EUVE. In this paper we re-examine 
the excess EUV emission in the Coma Cluster us- 
ing all the EUVE data available on this cluster. 
We obtain important new information on the char- 
acter of the EUV emission in this cluster. Given 
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these new results, we provide strong evidence that 
the EUV emission is produced by secondary elec- 
trons and positrons in the intracluster medium 
(ICM). This finding may well be the only secure 
evidence of the presence of secondary electrons 
and positrons in an intracluster medium. 

2. Data and Data Analysis 

All of the data employed were obtained with the 
Deep Survey (DS) telescope of EUVE (Bowyer & 
Malina 1991). In Table 1 wc provide an observing 
log of the observations. The total observing time 
was 390 ks. 

Various authors have used a number of data re- 
duction procedures in searches for EUV emission 
from clusters. Because of the misconceptions cre- 
ated by the use of incorrect analysis procedures, 
we describe the data reduction approach used here 
in some detail. These procedures were developed 
and documented in Bowyer et al. 1999. The valid- 
ity of this approach was examined and tested by 
Berghofer, Bowyer, & Korpela 2000 and its appro- 
priateness verified. 

First, the Coma data sets were screened to ex- 
clude noisy data. The pulse height distribution of 
each set was then examined and low energy counts 
produced by random noise were excluded by re- 
jecting counts below a low energy threshold. Since 
a low energy threshold is applied to the data by the 
onboard satellite data processing system, this step 
was not crucial. Indeed, Berghofer et al. 2000 have 
shown that changing the low energy threshold by 
as much factor of two has no effect on the end 
result. Nonetheless, this approach can improve the 
quality of the data set in at least some cases. Cos- 
mic rays interacting with the spacecraft and the 
detector produce a few high energy counts in the 
data which we removed by upper level threshold- 
ing. These counts are only a small fraction of the 



Date Duration (ks) 

12/25/95-12/28/95 50 
06/11/96-06/12/96 39 
01/12/99-01/14/99 53 
02/04/99-02/07/99 76 
03/15/99-03/21/99 172 

Table 1: Log of Observations 



total data set and ignoring this step does not sig- 
nificantly affect the end result. Nonetheless, these 
counts were easy to remove and we did so. Cor- 
rections were then made to account for telemetry 
limitations and detector dead time effects on the 
total observing time; these were ~ 10%. 

The next step in our analysis is quite important. 
A background was obtained from regions of the de- 
tector that do not view photons from the sky. This 
background arises from energetic charged parti- 
cles interacting with the satellite; these produce 
charged particles within the instrument that trig- 
ger counts in the detector. This background varies 
over time scales of weeks to months and depends 
upon geophysical conditions. Berghofer et al. 2000 
have shown that this background differs by only a 
factor of two over the course of the EUVE mission, 
but given the low counting rates from clusters of 
galaxies it is important that this background level 
be identified in order to establish the zero level for 
each particular observation. Accordingly, we es- 
tablished this background independently for each 
of the data sets wc employed. 

Most importantly, the correct telescope sensi- 
tivity profile, or exposure map, was used in con- 
nection with the analysis of the data. We note sim- 
ilar corrections for the instrument sensitivity over 
the field of view are routinely applied in the reduc- 
tion of most observations of diffuse X-ray emission. 
For example, observations of diffuse sources with 
the ROSAT PSPC are routinely corrected using an 
effective area exposure map (Snowden et al. 1994). 
Bowyer et al. 1999 have provided a map of the 
EUVE DS sensitivity profile using 363 ks of data 
from a variety of blank fields. The use of a sen- 
sitivity profile composed of a large number of in- 
dividual blank field data sets could, in principle, 
be questioned. Indeed, Lieu et al. 1999a claimed 
the EUVE DS sensitivity profile varies with time, 
but no analysis validating this claim was provided. 
Berghofer et al. 2000 carried out a detailed investi- 
gation of this possibility. They compared the 363 
ks data set referred to above with an assemblage of 
425 ks of data from a different set of blank fields 
obtained at different times. The two data sets 
were correlated at the 97% level, consistent with 
the statistical uncertainties in the counts in the 
individual cells in the two data sets. This demon- 
strated the stability of the EUVE DS telescope's 
sensitivity profile over time scales of years. In our 



2 



work on the Coma Cluster we used a sensitivity 
profile composed of 788 ks of data obtained by 
combining the two blank field data sets described 
above. 

Because of the different orientations of each of 
the different Coma observations, it was necessary 
to carry out the above steps on each of the individ- 
ual data sets separately. The results of each obser- 
vation were then summed. This required a knowl- 
edge of the absolute pointing of the spacecraft. Be- 
cause there are no obvious point sources that are 
present in all of the EUV images, it is non- trivial 
to confirm the pointing coordinates provided by 
the satellite. A comparison of the location of the 
maximum of the cluster emission in the images 
shows the relative pointing error in the nominal 
spacecraft pointing to be ~ 0'.28. Since this un- 
certainty is ^ the estimated point spread func- 
tion of the telescope, we simply added the images 
using the nominal spacecraft pointing. We note, 
however, that any conclusions based on the EUVE 
data will be uncertain at this, or smaller, scales. 

Next, the effects of absorption by the Galac- 
tic interstellar medium (ISM) on the EUV flux 
were determined. There are a number of programs 
available to determine the effects of the ISM on the 
X-ray flux from Galactic and extragalactic sources 
and any of these will provide a result that is es- 
sentially valid in the X-ray regime. However in 
the EUV, the situation is entirely different. In 
this band absorption is due only to hydrogen, neu- 
tral helium, and singly ionized helium. Metals can 
be ignored because they produce insignificant ab- 
sorption in comparison to these species, and the 
reduction of He I and He II due to the presence 
of He HI can be ignored because there is virtually 
no He III in the ISM (Heiles et al. 1996). The ap- 
propriate EUV cross sections must be used for H 
I, He I, and He II, and equally importantly, cor- 
rect columns are needed for each of these compo- 
nents. In particular, the amount of H II in the 
line of sight must be established in order to de- 
termine the true He I and He II columns. A full 
discussion of these issues and a comparison of the 
differing outcomes with the use of different com- 
pilations of cross-sections are provided in Bowyer 
et al. 1999. In this case we used a hydrogen column 
of 8.95x IQi'^ cm-2 (Dickey & Lockman 1990) with 
ionization fractions and cross sections for Galactic 
ISM absorption as described in detail in Bowyer 



et al. 1999. We note that Brcgman et al. 2003 have 
shown that small scale variations in the Galactic 
ISM can be as large as factor of 3 in some 1 degree 
fields containing clusters of galaxies and this can 
affect the magnitude of the EUV excess in these 
cases. However, this is not a factor in regard to 
the Coma Cluster where there is near spatial uni- 
formity of the Galactic Hi column as manifested 
in the NRAO map of this region with a spatial 
resolution of 21', and the finer scale IRAS 100 /xm 
map. 

We then derived the EUV emission produced 
by the high temperature X-ray emitting gas using 
Coma ROSAT PSPC archival data. We used a 
temperature of 9 keV (Briel et al. 1992). We note 
that a variety of temperatures, typically varying 
from 8 to 9 keV have been reported for the thermal 
gas in Coma by various authors. This variation 
has only a small effect on the ratio of the X-ray to 
EUV flux. This ratio for an 8 keV plasma is within 
10% of that for a 9 keV plasma. We used a factor 
of 128 to convert counts in the 0.5-2.4 keV band of 
the ROSAT PSPC to the EUVE DS-band counts. 
The PSPC conversion factor was derived from the 
MEKAL plasma code with abundances of 0.3 solar 
and a temperature of 9 keV. We corrected for the 
Galactic ISM as described above. 

The next task was to align the X-ray and EUV 
images. A source well away from the cluster cen- 
ter was detected at the same sky location in both 
the EUVE DS image and the ROSAT soft X-ray 
image. In both of these images the source was 
^ the point- spread functions of the respective 
detectors. A UV source, A 2305, is located within 
the central portion of the point-spread functions 
in both the EUVE and the ROSAT images. A 
QSO would typically produce a UV, EUV, and 
soft X-ray signature of this character. With this 
source as a fiducial, the images were aligned to 
^ 0'.43. We note that any comparisons between 
the EUV and X-ray data are uncertain at, or less 
than, scales of 0'.43. We then subtracted the EUV 
emission due to the X-ray gas from the total EUV 
emission detected by EUVE. 

The identification of diffuse emission in a sky 
map is difficult because of the low signal to noise 
ratio of data in individual cells of the map. This 
problem was identified early in the study of diffuse 
X-ray emission in clusters of galaxies. A solution 
univerally employed in studies of diffuse emission 
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Fig. 1. — (a) - The azimuthally averaged radial 
profile of the total 0.13 - 0.18 keV EUV count rate 
(solid), the EUV count rate in this band produced 
by the X-ray plasma plus the EUV background 
count rate (dashed), and the EUV background 
alone (dotted), (b) - The azimuthally averaged 
radial profile of the 0.13 - 0.18 keV EUV excess 
count rate in the Coma Cluster. 



in clusters is to construct the azimuthally averaged 
radial intensity profile of the flux. We derived this 
profile for the EUV emission in the Coma Clus- 
ter. The results are shown in Fig. 1. There is a 
substantial EUV excess out to ~ 14' and marginal 
evidence for emission to 20'. The dominant uncer- 
tainty in the determination of the overall excess is 
uncertainty in the determination of the ratio be- 
tween the EUVE and ROSAT c;ount rates due to 
the X-ray plasma. Including this uncertainty, the 
overall EUV excess is significant at greater than 
the 12(7 level. 

In order to obtain a value for the total EUV ex- 
cess in physical rather than instrumental units, we 



summed the excess counts shown in Figure 1 and 
computed the unabsorbed count rate by correcting 
for Galactic interstellar absorption as described 
above. We then divided by the EUV instrument 
effective area to obtain results in physical units. 
The unabsorbed EUV excess in the band from 68 
to 92 A (the approximate bandpass of the observa- 
tion as defined by the telescope high energy cutoff 
and the low energy cutoff produced by Galactic 
absorption) is 1.7 x 10""'^^ ergs~'^cm~^A~-'^ . As- 
suming a distance of 100 Mpc, this corresponds 
to a total energy output between 68 and 92 A of 
4.9 X lO^^ergs"'^. In this calculation we assmncd 
spectral indices between 1 and 1.6 which are ap- 
propriate given the source mechanism identified 
for the emission as discussed below. This result is 
relatively insensitive to the spectral index, with a 
variation in the flux of only a few percent for the 
index range listed. It is of interest to compare this 
energy output with the energy output of the X-ray 
plasma which is about 10^^ ergs/s, based upon a 
central density of 3 x 10""^ cm~'^ , a core radius of 
10'. 5, a /? of 0.75, a temperature of 9 keV, and the 
cooling function of Sutherland and Dopita 1993. 
If the EUV excess were due to a thermal plasma 
at 10^ K, the bolometric luminosity of this plasma 
would be 5 X 10^* erg/s, which is comparable to 
the energy output of the X-ray plasma. 

The value we obtain for the EUV excess, F\, 
is about a factor of two smaller than the number 
reported by Sarazin and Lieu 1998 after correct- 
ing for a difference in the assumed distance to the 
Coma Cluster. Although it is impossible to con- 
clusively identify the reason for this difference, we 
note that if we were to inappropriately compute 
the energy output using the full bandpass of this 
instrument rather than the cfFec'tivc bandpass, we 
would obtain a value similar to that reported by 
Sarazin and Lieu. 

3. The Relationship between the EUV Ex- 
cess and the X-ray Emission 

We first derived the ratio between the az- 
imuthally averaged EUV flux (0.13 - 0.18 keV) 
and the azimuthally averaged X-ray flux derived 
from the ROSAT 0.5-2.4 kcV X-ray band. We 
show this ratio as a function of increasing dis- 
tance from the cluster center in Figure 2. As can 
be seen this ratio is essentially flat. The error 
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bars increase at larger radii because of the limited 
EUV excess flux at these larger radii. 

An azimuthally averaged radial intensity profile 

is quite sensitive to the presence of diffuse emis- 
sion. However, by its very nature this process 
eliminates any possibility of examining details of 
the spatial distribution of the emission other than 
its average radial distribution. A direct study of 
the details of the EUV emission in this bandpass 
can only be achieved with a very substantial data 
set, which is now unobtainable. As an alterna- 
tive, we considered ways to investigate aspects of 
the spatial distribution which might prove to be 
useful. We first considered the number of EUV 
excess counts in individual cells in the sky map. 
The telemetered cell size of EUVE data is 4". 6. 
We summed these data into larger blocks. The 
minimum appropriate cell size is 0'.28 because the 
registration of the EUV images arc uncertain at 
this level. In addition, the use of a cell size smaller 
than the intrinsic resolution of the telescope could 
potentially provide misleading results. The re- 
sponse of the telescope is closely replicated by a 
Gaussian with a 90% included energy width of 1' 
and one possibility would be to convolve the data 
with a Gaussian of this size. However, we summed 
the counts in a 1 min square box since Hardcastle 
2000 has pointed out that the use of a smooth- 
ing function adds considerable uncertainty to the 



h ig. 3. Sky map oi the EUV excess m the Coma 
Cluster in 1' square bins (J2000). The isophotal 
lines shown are at 90%, 50% and 23% of the peak 
EUV excess emission. 



significance levels of the resultant data set. We 
then computed isophotcs of the EUV excess. The 
results are shown in Figure 3. The EUV excess 
appears to be more extended to the southeast al- 
though this is a region of low counts per bin. 

In Figure 4 we show the X-ray emission in 
the Coma Cluster derived from archival ROSAT 
PSPC data. A comparison of Figure 4 with Figure 
3 shows that the EUV emission is only detectable 
in the central portion of the X-ray image. This 
may be because the EUV emission is only present 
in the core, but it could simply be a sensitivity is- 
sue, and the lower intensity wings extend further 
out. We then carried out a standard linear cor- 
relation analysis between the EUV excess dataset 
shown in Figure 3 and the corresponding X-ray 
data shown in Figure 4. 

An immediate problem in carrying out a corre- 
lation analysis between the EUVE excess and the 
X-ray emission is that a correlation analysis will 
compare the number of counts in a given cell in one 
image with the number of counts in an identical 
cell in the other image without accounting for any 
statistical fiuctuations in these values. Hence with 
an image with a small number of counts in individ- 
ual cells, a false statement of a lack of correlation 
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Fig. 4. — X-ray emission in the Coma Cluster de- 
rived from ROSAT PSPC data. 



will be provided simply because of the statistical 
fluctuations of the data in the cells. 

To assess this effect for the EUV excess, wc car- 
ried out a correlation between two independent 
data sets of the EUV excess in the Coma Cluster 
as wc summed the coimts in individual cells into 
larger sized bins. The results are shown in Fig. 
5 as diamonds. As expected, the correlation be- 
tween the two independent data sets of the EUV 
excess is quite poor with smaller bin sizes, but in- 
creases rapidly as the bin sizes are increased and 
more counts are registered in each bin. 

We note that the determination of the confi- 
dence levels of the correlation measures shown in 
Fig. 5 arc inherently complicated by three statis- 
tical properties of correlation estimates: (a) they 
are inherently non-Gaussian, being mathemati- 
cally bomidcd to the interval -1 to +1; (b) they 
are asymmetrical; (c) their confidence intervals 
depend on the true population correlation value, 
which is unknown. This situation prompted Fisher 
1935 to create a nonlinear transformation of the 
correlation statistic into a Gaussian normal variate 
with uniform variance, namely Fisher's transfor- 
mation which depends only on the observed cor- 
relation value r and the number of independent 
data points A^, both of which arc known. In use, 
one simply converts the observed r into Fisher's z 
(which is Gaussian) and then converts the desired 
confidence interval in z back into an interval for r. 
Using this method we obtain the error bars shown 
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Fig. 5. — The correlation coefficients of two inde- 
pendent data sets of the EUV excess in Coma as 
a function of the size of the sky map binning is 
shown by diamonds. The correlation of the EUV 
excess with the X-ray emission as a function of the 
size of the sky map binning is shown by triangles. 



in Figure 5. We confirmed these error values by 

performing multiple Monte Carlo simulations of 
uncorrelated data. The standard deviation of the 
correlation coefficients of these uncorrelated sim- 
ulations in Fisher z-space was equivalent to the 
error values calculated using the number of inde- 
pendent points. 

We expected that there were sufficient counts in 
the deep ROSAT X-ray image that uncertainties 
in the photon statistics in the X-ray data would 
be inconsequential in comparison with the uncer- 
tainties due to the limited data in the EUVE data 
set. A self-correlation of the X-ray data verified 
this conclusion. 

We then carried out a correlation of the EUV 
excess image with the X-ray image as a function of 
increasing cell size. The results arc shown in Fig- 
ure 5 as triangles. Up to a scale of 4.0 arcmin^ the 
EUV self correlation and the EUV/ X-ray corre- 
lation both rise reflecting the limited quantity of 
EUV data. At larger scales the EUV self corre- 
lation is better than the EUV/X-ray correlation 
allowing us to make several definitive statements. 
First, there is a substantial, but not exact, spa- 
tial similarity between the EUV and the X-ray 
emission. This similarity peaks at scales of 8.8 
arcmin^ at a value of 0.86. There is an indica- 
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tion that the correlation falls off at larger scales, 
though this decrease is not significant at the three 
sigma level. We can also conclude that at scale 
sizes smaller than 4.0 arcmin^ the correlation is 
no better than 0.80 since a correlation cannot be 
improved by simply reducing larger image pixels 
to smaller pixels. 

In summary, we can state that at scales greater 
than 4.0 arcmin^ the spatial distributions of the 
EUV and the X-ray emission have substantial, but 
not exact, similarities. The correlation peaks at a 
value of 0.86 at a scale size of 8.8 arcmin^ with 
a suggestion that the correlation decreases at a 
larger scale. At scales less than 4.0 arcmin^ the 
correlation is no better than 0.80. 

4. Discussion 

In view of these new findings, it is useful to 
reconsider suggestions for the underlying source 
mechanism for the EUV excess in clusters of galax- 
ies. One proposal was that this emission was ther- 
mal emission from a "warm" (10^ K) gas (Lieu 
et al. 1996a; 1996b; 1999a,b; Mittaz et al. 1998; 
Bonamente et al. 2001). Buote 2000a, b carried 
out extensive analyses of ROSAT PSPC data and 
also claimed to have found evidence for gas at this 
temperature in the core of several clusters. 

The EUV emission in the Coma Cluster is al- 
most spherical and could conceivably be the prod- 
uct of a gravitationally bound gas. However, the 
maintenance of a warm intracluster gas is quite 
difficult to understand since gas at this tempera- 
ture is at the peak of its cooling curve and would 
typically cool in less than 0.5 Gyr (Landini & Mon- 
signori Fossi 1990). This has resulted in a fair level 
of skepticism in regard to a thermal origin for the 
EUV excess. A variety of observational studies 
have been carried out in an attempt to discover 
evidence for a warm 10^ K thermal gas. Initial 
studies with XMM- Newton showed no lines from 
a 10^ K gas in any of the clusters examined (Pe- 
terson ct al. 2001). However the Coma Cluster 
was not examined in these studies, which raised 
at least the possibility that Coma is uniquely dif- 
ferent and that the EUV excess in this cluster is 
indeed thermal in origin. Dixon et al. 2001 ob- 
tained long Far Ultraviolet Spectroscopic Explorer 
(FUSE) observations centered on the Coma clus- 
ter in search of Ovi 1032, 1038 emission which 



would be produced by a lO^K thermal gas. Al- 
though this emission was detected, the high res- 
olution of FUSE showed that all of this emission 
was Galactic and that none was red-shifted Ovi 
from material in the cluster. However, the pres- 
ence of thermal gas could not be completely ruled 
out since a large depletion of oxygen in the cluster 
would render these lines unobservable. 

Additional information on this topic has been 
obtained by Arnaud et al. 2001 and by Vikhlinin 
et al. 2001. Both of these groups studied the core 
of the Coma cluster in detail. Arnaud et al. ana- 
lyzed XMM-Newtori data and determined temper- 
atures in 3.5' X 3.5' regions in the 20' core of the 
cluster. The temperatures in these regions ranged 
from 7 to 8.5 keV (with a few outliers) with no 
evidence for lower temperature gas. Vikhlinin et 
al. used Chandra observations to search for lower 
temperature gas in very small regions in the core 
of the cluster. They found 1 to 2 keV gas within a 
7" radius of NGC 4874 and NGC 4889 which they 
attributed to emission from the halos of these in- 
dividual galaxies. However, immediately outside 
of these small regions they found only high tem- 
perature (9 keV) gas. 

Finoguenov et al. 2003 used XMM-Newton 
data in ~ 20' diameter bins and detected Ovii 
and Ovill emission ~ 30' off-center from the core 
of the cluster. These lines would be produced by 
a 0.2 keV (or 2 x lO'^ K) gas. They showed that 
this emission came from a filament in front of the 
Coma Cluster which was seen in projection against 
the cluster. The key difference between these mea- 
surements and the measurements of Vikhlinin et 
al. and Arnaud et al. was the fields of view in- 
volved. Only with the larger field of view is the line 
emission from the low temperature gas detected. 
Finoguenov et al. found that the oxygen line emis- 
sion was 1 / 30th of the X-ray emission of the 9 keV 
gas in the cluster center. Even if Finoguenov et 
al. were incorrect in their interpretation that this 
emission was from a filament in front of the clus- 
ter and was, instead, associated with 2 x 10^ K 
gas in the core of the cluster, its EUV emission 
would have been so faint as to be unobservable 
with EUVE and could not have been responsible 
for the excess reported here. 

Kaastra, Lieu et al. 2003 claimed to have found 
"warm" thermal emission within the central 12' of 
the Coma Cluster with XMM-Newton. However, 
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their claim is based on the (marginal) detection of 
a diffuse soft X-ray excess in the cluster and not on 
the detection of oxygen lines and it is contradicted 
by the work of Arnaud ct al. 2001 and Vihklinin 
et al. 2001. The Kaastra, Lieu et al. 2003 results 
could equally well be interpreted as non-thermal 
emission. 

Since the EUV flux is not thermal in origin, 
we have examined non-thermal processes as the 
source of this emission. Inverse Compton (IC) 
scattering of the 2.7 K cosmic microwave back- 
ground (CMB) photons by energetic electrons 
(Hwang 1997; Bowyer & Berghofer 1998; Enfilin 
& Biermann 1998; Sarazin & Lieu 1998) was sug- 
gested early on, and it is still the only suitable 
candidate non-thermal mechanism (Blasi & Co- 
lafrancesco 1999; Atoyan & Volk 2000; Brunetti 
et al. 2001a; Petrosian 2001; Tsay et al. 2002). 

Enfilin, Lieu & Biermann 1999 suggested clus- 
ter starlight radiation as the background photon 
field. The energy density in starlight photons is 
about two orders of magnitude less than that in 
the CMB in the core of Coma, so the starlight- 
IC model requires a number density of ~ 5 MeV 
electrons which is comparable to the thermal elec- 
tron density (Eni31in ct al. 1999). These particles 
then provide the dominant pressure in the cluster. 
This condition seems implausible both to estab- 
lish and to maintain. If such low energy electrons 
were mixed with the thermal plasma, they would 
transfer their energy on timcscalcs of a few hun- 
dred million years by way of Coulomb collisions 
which would result in excessive heating, even ig- 
noring heating by likely associated non-thermal 
protons. Magnetic fields strong enough to sepa- 
rate these non-thermal particles from the thermal 
plasma would lead to magnetic pressures exceed- 
ing the thermal gas pressure, which is similarly un- 
likely. If the posited 5 MeV electrons were relics of 
a much more energetic population, their original 
energy content would have been much greater than 
that already required for the starlight-IC model 
itself, making these problems worse. Thus, we re- 
ject this hypothesis as an luitenable explanation 
for the EUV excess in Coma. 

A number of authors have suggested specific 
IC-CMB models for the production of the EUV 
excess. Atoyan & Volk 2000 posited a popu- 
lation of "relic" electrons driven into the intra- 
cluster medium by galactic winds during intervals 



of galactic starbursts and then reaccelerated by 
strong merger shocks. To avoid excessive radio 
emissions from the assumed hard energy spectrum 
of the EUV electrons (sg = 2.1) with multi-/iGauss 
magnetic fields, their model included an ad hoc 
electron energy cutoff near 250 MeV. Brunetti 
et al. 2001b proposed a model for EUV IC- 
CMB emissions based on turbulent reacceleration 
of non-thermal electrons recently injected by the 
head-tail radio galaxy NGC 4869, which is sev- 
eral arcminutes west of the cluster center. This 
was an extension of a model that would explain 
the radio emission and the 40 KeV X-ray emis- 
sion in the Coma Cluster (Brunetti et al. 2001a). 
Important constraints imposed were that this pop- 
ulation should not produce (an unobserved) 7-ray 
emission, and it should produce the observed spec- 
tral steepening of the radio emission with increas- 
ing distance from the cluster core. This model in- 
cluded an initial cosmic ray population produced 
~ 2Gyr ago, reacceleration ~ IGyr ago by shocks 
from mergers, and a recent injection of low energy 
cosmic rays that is responsible for the EUV emis- 
sion. Their model assumed a relatively hard elec- 
tron spectrum for the EUV electrons (hereafter, 
EUVe) with Sg = 2.6, and a magnetic field ~ 0.5 
/iGauss, and required a cutoff near 500 MeV to 
avoid emissions conflicting with observations in 
other bands. The relative complexity of the full 
model illustrates the difficulty in finding a unified 
model for non-thermal emissions in Coma. 

Sarazin & Lieu 1998 proposed a model for EUV 
emission in clusters in which relic low energy elec- 
trons accumulated from various origins would be 
distributed similarly to the thermal plasma, that 
is HEuve oc Ute- That distribution predicted an 
azimuthally averaged ratio, Ieuv/Ix l/nte, in- 
creasing substantially with distance outside the 
cluster core. Berghofer et al. 2000 derived this 
ratio for the Virgo cluster as a test of this model. 
They found that the ratio was flat with increasing 
distance from the center of the cluster in contra- 
diction to the prediction of this model. We have 
derived this ratio for the Coma Cluster, and, as 
shown in Figure 2, it is essentially constant. Us- 
ing a beta model for Coma (Briel et al. 1992) with 
/3 = 0.75, Tc = 10'. 5, and assuming nsuve oc rite 
one finds that Ieuv /Ix should have increased by 
about a factor of two from the cluster center to a 
distance of 10'. 5 and by a factor of six at 21'. All 
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of these outcomes are clearly inconsistent with the 
data. 

Three previous studies considered secondary 
emission in connection with the EUV excess. Blasi 
& Colafrancesco 1999 considered secondary emis- 
sion as part of a unified model for non-thermal 
emissions in Coma. They found their model had 
multiple problems. The spatial distribution of 
the radio emission was not correct, and too much 
gamma radiation was produced. Finally, the EUV 
emission produced was too low. Blasi 2001 mod- 
eled emission from secondaries as part of a treat- 
ment of nonthermal emission in cluster mergers, 
including Coma. He assumed that a strong merger 
shock would inject primary electrons and protons 
with a density distribution proportional to the 
thermal plasma. With the parameters he em- 
ployed, the associated EUV emission would be 
dominated by primary electrons, so the resultant 
EUV spatial distribution would take the same 
form as that proposed by Sarazin and Lieu (1998), 
and would be incompatible with the observational 
results reported here. Miniati et al. 2001b esti- 
mated the EUV flux from secondary emission in 
clusters as part of a larger study of cluster forma- 
tion. Their EUV luminosity vs cluster tempera- 
ture relations underestimated the observed EUV 
flux in Coma by about an order of magnitude. 
Consequently, Miniati et al. did not pursue the 
idea that IC-CMB emission from secondary elec- 
trons could be the underlying source mechanism 
for the EUV excess. In retrospect, their luminos- 
ity estimates were artiflcially low because of the 
effects of flnite numerical resolution in their simu- 
lations. This signiflcantly reduced the central gas 
densities in clusters which resulted in an underes- 
timate of the secondary emission flux. 

No existing models produce an EUV inten- 
sity distribution that is highly correlated with the 
thermal X-rays as shown in Figure 5, and simul- 
taneously produce a constant ratio between the 
azimuthally averaged EUV and X-ray intensities 
as illustrated in Figure 2. Therefore, we have 
searched for a new model that would naturally 
produce these outcomes. 

4.1. General Constraints 

Before introducing a specific model that will 
yield these observational findings, we first estab- 
lish some general constraints on the emitting par- 



ticles and their environment that would apply 
to any successful model. In the IC-CMB sce- 
nario, the EUV excess is produced by electrons 
of characteristic energy E ^ 200^ei50ey MeV 
(7 ~ 400^ei50ey), where eisoeV is the EUV pho- 
ton energy, normalized to 150eV (A ~ 80 A). 

The magnitude; and distribution of the mag- 
netic field in the cluster are important constraints 
in any model for non-thermal emissions in Coma. 
Extensive work has been carried out in efforts 
to determine the magnetic field strength. Re- 
cent summaries of the observational situation and 
possibility of reconciling the (apparently) contra- 
dictory results have been provided by Kronberg 
2003; Clarke 2003; and Brunetti 2003. Differ- 
ent approaches yield different results. One ap- 
proach is to calculate the field based on the as- 
sumption of equipartition between the energy den- 
sity of the relativistic particles associated with the 
radio emission and the magnetic field. A recent 
and especially detailed result for the Coma Cluster 
using this approach has been obtained by Thier- 
bach et al. 2003 who find an equipartition field of 
~ 0.7/xGauss if electrons are the relativistic gas, 
or ~ 1.9/xGauss if the proton-to- electron energy 
density ratio in the relativistic gas is the same as 
that in the ISM. Faraday rotation measures of ra- 
dio sources in clusters have been extensively stud- 
ied as a means of determining cluster magnetic 
fields. Very high fields have been obtained us- 
ing sources embedded in the cores of clusters with 
cooling centers (EUek 1999; Taylor et al. 1999). 
Rotation measures of radio sources behind clus- 
ters have been measured by a number of groups 
(Kim et al. 1990; Feretti et al. 1995). The most 
extensive results using this approach have been 
obtained by Clarke et al. 2001. They find fields 
that are typically in the range of 5 to 10 /xGauss. 

Fields in the range of 0.1 to 5 /iGauss are re- 
quired if the EUV excess is the product of IC- 
CMB in a uniform magnetic field (e.g. Hwang 
1997; Atoyan & Volk 2000; Brunetti et al. 2001b; 
our discussion below), but the higher values in 
this range can only be realized if a rather arbi- 
trary high energy cutoff is imposed upon the un- 
derlying cosmic ray spectrum, or if the cosmic 
ray spectrum is very steep. A less extreme ex- 
planation for the lower fields required in IC-CMB 
models for the EUV is that the magnetic fields 
are not homogeneous and that the EUV excess 
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originates in low-field regions while high-field re- 
gions produce the higher Faraday rotation mea- 
sures (Petrosian 2001; Newman et al. 2002; Beck 
et al. 2002). Trcgillis ct al. 2003 studied synthetic 
radio and X-ray images derived from high reso- 
lution three-dimensional MHD radio galaxy sim- 
ulations to compare average field estimates with 
actual magnetic field properties in the simulated 
objects. They found that the estimated fields 
roughly corresponded to actual rms fields, but 
scattered around the physical rms value by a fac- 
tor of ~2-3. In light of the above discussion, we 
assume a field of about 1 /iGauss in the ICM of 
the Coma Cluster in the following. 

Energy loss timescales provide an important 
general constraint on models for the EUVe. For 
conditions in the X-ray core of Coma, approxi- 
mately coinciding with the EUV excess (thermal 
electron density, rite « 3 x 10~^cm~^, Briel et al. 
1992), IC-CMB and Coulomb energy losses are 
currently roughly comparable at 200 MeV (e.g. 
Sarazin 1999; Petrosian 2001). However, since 
synchrotron energy losses compare to IC losses by 
the ratio {B/B^}'^ and cx (1 + z)^, this ratio 
was probably smaller in the past. With a magnetic 
field of about 1 /uGauss, synchrotron losses can be 
neglected. Using standard relations (e.g. Sarazin 
1999: Petrosian 2001) and correcting the IC-CMB 
loss rate for the Hubble expansion with go = 0.5, 
but assuming rite was not greater in the past, it 
is simple to demonstrate that IC-CMB losses were 
dominant over Coulomb losses at these energies, 
and that the characteristic IC-CMB energy-loss 
lifetime for EUVe is ~ 2 Gyr. 

Another key point is that electrons at these 
energies diffuse very slowly in the intracluster 
medium (Schlickeiser et al. 1987; Volk et al. 1996) 
For Bohm diffusion in a /iGauss field, 200 MeV 
electrons would diffuse only about 10 pc during 
their lifetimes. In almost any plausible cluster 
field and turbulence model the EUVe are effec- 
tively tied to the local plasma. Mixing of the 
cluster plasma will take place on timescales of Gi- 
gayears in response to stirring in the cluster (e.g. 
Markevitch, Vikhlinin & Mazzotta 2001) caused 
by mergers and AGN activity. But since the 
EUVe electrons are tied to the cluster medium, 
the spatial distribution of both of these species 
will be similar, though not identical. In clus- 
ters with current active energy deposition such 



as Virgo or Hydra, fresh particle populations will 
not become mixed immediately, as illustrated by 
the X-ray holes seen in such clusters (e.g. Nulsen 
ct al. 2002). On Gigayear timescales, however, the 
nonthermal particles will become mixed through- 
out the cluster. 

We next explore spectral constraints on the 
EUVe population that can be derived by requiring 
that it does not produce emission in other bands 
that exceeds those observed. The EUVe popula- 
tion directly includes only energies near 200 MeV, 
but is likely to continue to higher energies follow- 
ing a normal power-law spectrum. In particular, 
we see no reason to introduce an artificial cutoff at 
higher energies. Clearly an important constraint 
is that the high energy extension of the EUVe 
population does not produce IC-CMB in excess 
of the observed non-thermal high energy excess 
(hereafter HRX) . Both BeppoSax (Fusco-Femiano 
et al. 1999) and RXTE (Rephaeli & Gruber 2002) 
gave results that can be expressed in terms of a 
flux near 40 keV of z/F,y ss 8 x 10"-^^ erg cm~^s~-^ 
inside a radius ^ 1 degree. If the EUVe pop- 
ulation includes electron energies approaching 4 
GeV, its IC-CMB spectrum will reach into this 
band. Our measured EUV flux corresponds to 
uF^ = 1.4 X 10~^^ erg cm~^s~^ at 150 eV. A sim- 
ple power-law extension of the IC-CMB spectrum 
from 150 eV to 40 keV with a spectral index, a, 
{Fv oc corresponding to an electron energy 

distribution ne{Ee) oc E~^^, where Se = 2a + 1) 
would fall below the observed excess 40 keV fltix 
if a > 1.1, or Se > 3.2. Since the HRX field is 
substantially larger than the EUV source, the ex- 
tended spectrum could be steeper, so Sg > 3.2 is a 
conservative limit above Eg ~ 200 MeV. 

Similarly, extension of the EUVe population to 
higher energies could contribute detectable radio 
emission. This constrains both the form of the 
electron spectrum and the effective magnetic field 
strength, as discussed previously. The lowest fre- 
quency synchrotron flux measured for Coma C is 
49 Jy at 30.9 MHz (Giovannini et al. 1993). The 
radiating electrons would have characteristic ener- 
gies Ee ~ 1.4 GeV B~^^^, where B± = B{cose) is 
the sky-plane component of the source magnetic 
field expressed in /xGauss. For B± ^ 1 /iGauss, 
these electrons would be roughly an order of mag- 
nitude more energetic than those producing the 
EUV emission. The requirement that the observed 
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radio flux exceed any synchrotron flux, F,y^ , pro- 
duced by a high energy extension of the EUVe 
population can be conveniently expressed by the 
constraint Rst = VsP^JviFy^ = \sFxJ\F-x. < 
1.07 X 10-3, where F^, is the observed IC-CMB 
flux in the EUV. 

Assuming a power-law electron energy distri- 
bution over the relevant range, and that the EUV 
emission is IC-CMB, the ratio of the associated 
radio synchrotron flux to the EUV flux in a uni- 
form magnetic field is easily shown to be (Jones 
et al. 1974) 



B± = 



■be / \ "^1 



3 jao \Vil'fiB 



Rs 



(1) 



where Vg = 30.9 MHz and Ui = 37 PHz are the 
observed frequencies for synchrotron and IC-CMB 

emission, jao ~ 1 and ~ 1 are constants tab- 
ulated in Jones et al. 1974, = kTcMs/h = 
57 (1 + z) GHz, = 3.2 {l + zf ^Gauss is a fidu- 
cial field strength whose energy density matches 
the CMB, and i/^s = eB^/{2-Kmc) = 9 (H- zf 
Hz. 

A range of values for the integrated Coma C ra- 
dio spectrum is viable; the appropriate choice de- 
pends on the (unknown) underlying source model 
(Thierbach et al. 2003). Reasonable values for the 
spectral index range from a = 1.35 to 0.83. For 
consistency, values of Rsi near the observed limit 
would require model spectra at least as steep. The 
associated magnetic field constraint from equation 
1 depends on the spectrum chosen. For example, 
it would range from B±^ ~ 0.1 /xGauss for a = 1.0 
{sp = 3) to B± ^ 1.4 /iGauss for a = 1.5 (sg = 4). 
If the magnetic field is isotropically oriented, these 
values of B± should be increased by roughly 20% 
to arrive at an estimate of the magnetic field. 

4.2. Secondary Electrons as the Source of 
the EUV Flux in Coma 

We now propose a specific model that pro- 
duces the observational results presented here, and 
then demonstrate that this model does not vio- 
late the general constraints derived above. The 
observed correspondences between the EUV and 
thermal X-rays strongly suggest that the low en- 
ergy EUV emitting particles and the thermal intr- 
acluster medium have a physical interdependence. 
As pointed out by Sarazin and Lieu 1998, the 



intuitively obvious relation between cosmic rays 
and thermal plasma is one in which the cosmic 
ray particles are relatively well mixed with the 
thermal plasma, i.e. in which Uc cx n^. This 
would be a natural expectation if the cosmic rays 
were accumulated over much of the formation of 
the cluster, and its various constituent compo- 
nents mixed following mergers. A scenario that 
incorporates this expectation and leads to the re- 
quired observational outcome of nsuve nf^ is 
that the EUV emission is the result of secondary 
electrons and positrons (hereafter SEP) produced 
as byproducts of inelastic collisions between well- 
mixed primary cosmic ray protons and the ther- 
mal intracluster plasma. The production rate for 
SEP scales as ricp X ntp, where ricp is the den- 
sity of cosmic ray protons. Then if the SEP en- 
ergy losses are independent of cluster position, as 
they would be for IC-CMB-dominatcd losses, the 
consequent distribution of EUVe is the required 
riEUVe ^ nsEP oc nfp nfg. 

A variety of primary cosmic ray source mech- 
anisms such as supernovae, active galaxies, and 
terminal galactic wind shocks might lead to the 
needed spatial distribution for this explanation of 
the EUV excess, so long as their contributions 
were spread over enough time and were spatially 
distributed into the various components. 

However, another source seems to us to be 
the most likely candidate for these particles: cos- 
mic rays accelerated at the large scale "structure 
shocks" that accompany cluster formation. These 
include what are usually termed merger and ac- 
cretion shocks, although recent cosmology simu- 
lations demonstrate a more complex and richer 
shock pattern than those labels suggest (Miniati 
et al. 2000; Ryu et al. 2003). In a cosmic structure 
simulation that included shock-accelerated cosmic 
rays Miniati et al. 2001a indeed found that the cos- 
mic ray proton distribution in the central regions 
of their clusters scaled roughly with the thermal 
gas although the two distributions did show dif- 
ferences from cosmic rays associated with recent 
shocks, especially outside the cluster core region. 

A more detailed evaluation of our model re- 
quires that we estimate the population of cosmic 
rays responsible for the SEP. Inelastic collisions 
between cosmic ray protons with kinetic energies 
above about 300 MeV and the thermal intraclus- 
ter medium produce mainly charged and neutral 



11 



pions. The charged pions decay into muons and 
neutrinos, and the muons into the SEP that are 
responsible for the EUV emission. We express the 
cosmic ray proton density distribution as a power- 
law of the form 



n(7p) 



(2) 



where Ep = ^ypMpC? is the proton energy. 

The approximate formalism given in Mannheim 
and Schlickeiser 1994 gives us a simple expression 
for the net omnidirectional SEP production rate, 
^e, with energies jmeC^ >> 35 MeV from p-p col- 
lisions. In particular, for a power law proton en- 
ergy distribution we have 
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/ Mp 



Qe « Y^appCntpHpo 



7 cm ^s \ 



(3) 

where (Tpp Ri 3.2 X 10 ^^cm ^, Sgo = ^{sp — ^), a,nd 
it 9e ■ 1*^^ energies this expression overes- 
timates SEP production, so we have compared it 
with a numerical calculation for qf based on the 
more acciirate pion production described in asso- 
ciation with equation 5 (following) and the SEP 
distribution given by Moskalenko & Strong 1998. 
For electron energies of 150 MeV the two results 
agree to better than about 50% so in what follows 
we use the simpler expression in equation 3. We 
argued above that the lifetimes of the EUVe at 
~ 200 MeV SEP are determined by their energy 
losses against IC-CMB. We can reasonably assume 
that the SEP density is set by a balance between p- 
p production and IC-CMB losses. Since the EUV 
emission is the same IC-CMB, it is straightforward 
to derive the expected omnidirectional EUV vol- 
ume emissivity directly in terms of the cosmic ray 
proton density. The result using equation 3 is 



Jcto 



26 



12(sp-5/4) \2Am, 



(4) 



appCntpUpoirieC [ — 

^BC 



erg cm ^ s ^, 



where Vi , f ^ and j^^ were identified in relation to 
equation 1, Sg = |(sp + 3) = Seo + 1 is the spectral 
index of the steady-state SEP energy distibution, 
and, once again, a = (sg — l)/2. With a con- 
stant cosmic ray density fraction, fp = Upo/ritp, 
the EUV emissivity scales with nj , as required 



Using this result we can integrate over the clus- 
ter to compute an EUV luminosity, which can then 
be compared with the observational result. The 
result will depend on an assumed cosmic ray en- 
ergy spectral index as well as on the intracluster 
medium density distribution. If the cosmic ray 
proton flux is due to structure formation shocks, 
the cosmic ray spectrum represents an average 
from the shocks dissipated in the local gas over 
cosmic time. In the test particle limit for diffusive 
shock acceleration, which is a reasonable approx- 
imation for relatively weak shocks, the standard 
relation is Sp = 2{M'^ + 1)/{M'^ - 1), where M is 
the shock Mach number. Strictly speaking, this 
index applies to the momentum spectrum, n{pp), 
of cosmic ray protons accelerated at shocks, rather 
than the energy spectrum, n(7p), that we defined 
in equation 2 in order to apply analytic expres- 
sions for our simple model estimates. The two 
spectra compare as 'n{'jp) / n{pp) = (1 — 1/7^)*'"/^. 
At relativistic energies the two forms converge; 
at the threshold for pion production, 7p w 1.3, 
n(7p)/n(pp) « 1/3, which roughly compensates 
for the overestimate in SEP production from equa- 
tion 3. 

In a recent detailed analysis of shocks formed in 
a high resolution cosmic structure formation simu- 
lation, Ryu et al. 2003 found that the most impor- 
tant shocks for cosmic ray acceleration were those 
with M 2 — 4, corresponding to Sp ~ 2.3 — 3. 
To be specific in our estimates below we choose 
Sp = 2.5, which leads to Sg = 11/3 and a = 4/3 
which will be consistent with the radio halo and 
HXR constraints. Then assuming a beta law dis- 
tribTition for the intracluster medium, ntp{r) = 
nio/(l + (r/a)2)3/3/2, ^ith rito = 3 X lO'^ cm'^, 
/3 = 0.75, (Briel et al. 1992) and a = 300 kpc 
(corresponding to 10'. 5 at 100 Mpc), we compute 
the spatially integrated EUV flux to be UiF^. w 



1 X 10 /„ erg cm 



Comparing this result to 



by the EUVE data for Coma. 



the observed I'iF^. 1.4 x 10 erg cm ^s ^, we 
obtain /p « 1.4 x 10"''. 

A test of the reasonableness of this result can 
be made by determining its consistency with up- 
per limits on the 7-ray flux in the Coma cluster. In 
addition to charged pions, inelastic p-p collisions 
will produce neutral pions, which will quickly de- 
cay to 7-rays. The resultant 7-ray spectrum peaks 
near 70 MeV, but extends to higher energies and, 
in particular, into the > 100 MeV EGRET band. 
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The 7-ray cmissivity due to a power-law cos- 
mic ray spectrum approaches a power-law at high 
energies, making it relatively straightforward to 
compiitc an analytical estimate of the high energy 
flux. However, the EGRET band is too close to 
the 70 MeV peak for that approximation to be ad- 
equate for our needs. Fortunately, semi-empirical 
relations for our range of interest are available in 
the literature. We have followed the formulation 
laid out conveniently in Schlickeiser 2002. The 
omnidirectional 7-ray emissivity can be written as 



(5) 

where q„o{E„) ex CTppCntpUpo is the neutral pion 
production rate. That rate asymptotes to a power- 
law at high energies with spectral index ,s,.o, jnst as 
for the SEP, but drops sharply as the pions become 
nonrelativistic. Using the full expressions given 
by Schlickeiser, assuming Sp = 2.5 and integrating 
equation 5 over photon energy, we obtain [E^ > 
lOOMeV) K, 0.32appCntenpo photons cm~^ s~^. 
From the values established above for Coma, in- 
cluding the required cosmic ray density fraction, 
fp, we obtain an estimated 7-ray flux Fj{Ej > 
lOOMeV) w 1.4 X 10"*' photons cm~^ s"!. 
Sreekumar et al. 1996 give a 2 a upper flux limit 
in this band of 4 x 10~^ photons cm~^ s~^ for the 
Coma cluster. This is well above the 7-ray flux 
produced by our required SEP population. 

Simulations such as those of Miniati et al. 2000, 
2001a, b and Ryu et al. 2003 have suggested that 
structure formation shocks might lead to cosmic 
ray energy pressures approaching as much as 1/3 
the total intracluster nicidium pressure. Accord- 
ingly, we have estimated the cosmic ray energy 
density in the core of Coma that would be neces- 
sary under the SEP model we have proposed. For 
the proton spectra of immediate interest {Sp ~ 
2.5), most of the kinetic energy resides in mildly 
relativistic particles, independent of whether we 
use the energy power law of equation 2 or the 
analogous momentum power law. We can write 
%poMpc'^, which gives Ucp ^ 



approximately Uc^ 

10~^^erg cm~^ for the nonrelativistic cosmic ray 
proton distribution. This compares to the ther- 
mal energy density, SritpkT w 4 x 10~'^^erg cm~^. 
This rough estimate can be compared to estimates 
from cosmology simulations which range upwards 
of tens of percent (e.g. Miniati et al. 2001a; Ryu 



et al. 2003). 

We can summarize the constraints on our model 
as follows. The electron energy spectrum must 
have a power law slope steeper than Se ~ 3-2, in 
order to avoid excess inverse-Compton hard X-ray 
emission. Constraints set by the observed radio 
synchrotron flux depend on the spectral index of 
the radio emission and the cluster magnetic field. 
An electron spectrum with Sg = 3, requires the 
cluster magnetic field be ~ 0.1/xGauss. If Se = 
4 the cluster field must be ~ 1.4/iGauss. These 
outcomes are shown graphically in Figure 6. 

5. Conclusions 

We have analyzed archival data obtained with 
EUVE on the core of the Coma cluster. We find 
the ratio of the azimuthally averaged EUV and 
X-ray intensities is essentially constant with in- 
creasing cluster radius. In addition, a correlation 
of the diffuse EUV emission with the diffuse X-ray 
emission shows that the detailed spatial distribu- 
tions of these emissions are quite similar, but not 
identical. 

XMM-Newton and Chandra observations show 
that there is no intracluster 10^ K gas in the 
core of the cluster and hence that the EUV emis- 
sion must be non-thermal. The only viable non- 
thermal source for the EUV excess is photons in- 
verse Compton scattered by 200 MeV electrons 
from the cosmic microwave background. To ac- 
count for the observed EUV intensity distribu- 
tion, the scattering electrons must be distributed 
with a spatial density roughly in proportion to the 
square of the thermal plasma density. A scenario 
that naturally produces these results is that the 
EUV emitting particles are electron/positron sec- 
ondaries produced by inelastic collisions between 
primary cosmic ray protons and thermal protons 
with similar spatial distributions. This type of dis- 
tribution would be a natural outcome if the cosmic 
rays had been produced over the past several Gyr 
and had become well-mixed with the intracluster 
medium. 

This model accounts naturally for the average 
spatial distribution of the EUV emission in rela- 
tion to the X-ray emission. It accounts for the sim- 
ilar, but imperfect, pixel-to-pixel correspondence 
of the EUV and X-ray emission. Finally, it accom- 
plishes this without violating observational limits 
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log(hi/ (eV)) 



Fig. 6. — A plot showing observations constrain- 
ing the secondary emission model. The solid lines 
are a schematic representation of the radio obser- 
vations with the range of spectral indices justi- 
fied by Thierbach et al. 2003. The solid squares 
are observed values of the EUV and HRX fluxes. 
The open square is the observational upper limit 
to the gamma-ray flux. Our model produces the 
observed EUV flux while not exceeding the radio 
and HRX fluxes for a range of reasonable mag- 
netic fields. The open triangle shows the 7-ray 
flux produced by our model. 



in other bands of the spectrum. We have demon- 
strated that the required underlying cosmic rays 
could reasonably have been produced in large scale 
structure shocks accompanying the cluster forma- 
tion. 

Secondary electrons as the source of the emis- 
sion in radio halos were first suggested by Denni- 
son 1980 and have since been discussed by a large 
number of authors. However, observational evi- 
dence for the presence of these particles has been 
lacking. There is, in fact, a growing body of ev- 
idence that, in general, these particles may not 
be the underlying source of the emission in radio 
halos (Brunetti 2003; Kuo et al. 2003). Hence 
the EUV emission in the Coma cluster may be the 
only direct evidence for secondary electrons in an 
intracluster medium. 
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Note added in proof.- Miniati has recently remodeled the nonthermal emission from clusters with an 
improved treatment (F. Miniati, MNRAS, 342, 1009 [2003]). In his published work, he shows results only 
above 10 keV. However, he has now extended his work to lower energies and has calculated the EUV flux 
with the specific parameters of the Coma Cluster (F. Miniati, 2004, private communication). He finds that 
the IC EUV emission from SEPs is consistent at the 20% nominal level with the measured EUV flux reported 
here. Further the spatial distribution for this emission is concentrated in the inner Mpc of the cluster. These 
results provide additional support for the ideas presented here. 
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